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Abstract 

We calculate the transverse momentum pT distribution for production of the T states in hadronic 
reactions. For small pt{^ My), we resum to all orders in the strong coupling the process- 
independent large logarithmic contributions that arise from initial-state gluon showers. We demon- 
strate that the pT distribution at low pj- is dominated by the region of small impact parameter b 
and that it may be computed reliably in perturbation theory. We express the cross section at large 
Pt by the 0{al) lowest-order non- vanishing perturbative contribution. Our results are consistent 
with data from the Fermilab Tevatron collider. 
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I. INTRODUCTION 



The theoretical description of the transverse momentum pt distribution of heavy quarko- 
nium production in hadron colhsions raises interesting challenges. Most calculations within 
the framework of perturbative quantum chromodynamics (QCD) consider the distribution 
at large pt at collider energies and tend not to address the region of low pt where the cross 
section is greatest and the bulk of the data lie . A purely phenomenological fit to 

the low Pt data on T production |J| appears to require sizable non-perturbative parton- 
krp smearing From a theoretical point of view, the region of low px is expected to be 
influenced strongly by initial-state gluon showering. A flxed-order perturbative treatment 
in QCD leads to singular terms in the region of small px of the type l/p^, enhanced by 
large higher-order logarithmic contributions caused by initial-state gluon radiation. These 
contributions have the form aslog^{My/p'^) for every power of the strong coupling a^, and 
reliable predictions, especially in the regions of small and moderate pt, require that the 
logarithmic contributions be summed to all orders in a.^- 

The Collins-Soper-Sterman (CSS) impact parameter 6-space resummation formalism 
has been used successfully for the all-orders resummation of large initial-state logarithmic 
terms in several cases of physical interest 
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14| . In this paper, we argue 



and demonstrate that the resummation formalism should apply at the scale of the T mass 
in hadronic collisions at collider energies. We extend the formalism and use it to compute 
the Pt distribution of the T states jisl . We obtain good agreement with the data jj, on 
the Pt distribution of T production at Tevatron energies for all pt- 

Different from the production of the W, Z, and Higgs bosons, or of a virtual photon 
in the Drell-Yan process, the T is unlikely to be produced in pointlike fashion in a short- 
distance hard collision. Instead, a bottom quark bb pair is produced in the hard collision 
and then transmutes into a colorless T meson through soft radiation and coherent self- 
interaction. Therefore, there are issues to address before the CSS formalism can be applied 
to T production. These include the color structure of the lowest order production processes: 
qq bb{Q) and gg bb{Q), and the relatively small value of the bb pair mass Q. 

Most applications of the resummation formalism are to the production of systems that 
are color singlets whereas the bb system produced in qq bb and gg —>■ bb need not be 
color neutral. Nevertheless, because the b quark mass is large, gluon radiation is suppressed 
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from the final-state heavy quark hnes and from virtual exchange lines that lead to the pro- 
duction of heavy quark pairs Correspondingly, the important logarithmic terms are 
associated with gluon radiation from the active initial-state partons, the same as those in 
massive lepton-pair (Drell-Yan) and Higgs boson production. The process-independent lead- 
ing logarithmic terms do not depend on the color of the heavy quark pair. Color dependence 
becomes relevant for the higher order terms, as explained in Sec. IIVI 

The overall center-of-mass energy dependence of the CSS 6-space distribution function 
is examined by Qiu and Zhang in Ref. j^. They show that the location of the saddle point 
of this distribution can be well within the perturbative region of small b for Q as small 
as 6 GeV at the Tevatron collider energy. The resummed 6-space distribution is peaked 
strongly in the perturbative region of small b, as we show in Sec. IIVI and the px distribution 
of T production should be amenable to a resummation treatment. Despite the fact that the 
logarithmic term In Q is not large, the large value of opens a large region of phase space 
for gluon emission. Correspondingly, as is demonstrated in this paper, the shape of the pt 
distribution for T production is determined by the resummable part of the gluon shower 
and is predictable quantitatively at low px- 

We begin in Sec. |n] with the basic assumption that the px distribution of T production 
is derived from the px distribution for the production of a pair bb of bottom quarks. We 
express the differential cross section in terms of a two-step factorization procedure. We 
present our fixed-order perturbative calculation applicable at large transverse momentum 
in Sec. IIIII where we also describe models that specify the manner in which the bb pair 
transforms into the T. In Sec. II V| we specialize to the situation at small px and summarize 
the required parts of the all-orders resummation formalism. Section |3 is devoted to our 
numerical results and comparison with data. We provide of our conclusions and discuss 
potential improvements of our calculation in Sec. IVIl 

II. PRODUCTION DYNAMICS 

We use a two-step factorization procedure to represent production of the T states, with 
particular attention to the prediction of transverse momentum distributions. We begin with 
the assumption that a pair of bottom quarks bb is produced in a hard-scattering short- 
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distance process: 

A{pa) + B{pb) bb{Q) H T(p) +X]+X'. (1) 

Because the mass Q of the bb pair is large, the pair is produced at a distance scale ~ 
1 / {2mf,) ~ 1/45 fm. This scale is much smaller than the physical size of a T meson. The 
compact bb pair may represent the minimal Fock state of the T, but the overlap of this 
minimal Fock state with the full wave function of the T is perhaps small, as is suggested 
by the inadequacies of the color-singlet approach 1^1 in some situations 1^1, and other 
components of the wave-function must be considered. Alternatively, one may realize that 
the compact bb system is unlikely to become an T meson at the production point. Instead, 
the pair must expand, and the b and b will interact with each other coherently until they 
transmute into a physical T meson. 

Once produced in the hard-scattering, a bb pair of invariant mass Q > 2Mb is more 
likely to become a pair of B mesons. Therefore, the virtuality of the intermediate bb pair 
should be limited if an T is to result. This limitation of the virtuality allows us to use 
perturbative factorization and to write the differential cross section in the usual way Q, Q| . 
For pt ^ 2{Mb — frif,), we write 



^^^^ = T.jd^a <Pa/A{x.) dx, ^^^^ . (2) 



In Eq. Pt and y are the transverse momentum and rapidity of the final T. The functions 
(t)i{x) are parton distribution functions; Xa and x^ are fractional light-cone momenta carried 
by the incident partons; and Eq. Q expresses initial-state collinear factorization. The 
spectator interactions between the beam remnants and the formation of the T meson are 
suppressed by one or more powers of l/p^. 

Since the momentum of heavy quark b (6) in the pair's rest frame is much less than the 
mass of the pair, Q — 2mh < 2Mb — 2mf, <^ 2mb, the parton-level production cross section 
ddah^rxl dp\dy in Eq. ^ might be factored further as is sketched in Fig. [T] The 
incident partons labeled Xa and Xh interact inclusively to produce an off-shell hb system plus 
state X' . In turn, the hb system evolves into the T plus a system labeled X\ X = X' + X 
This second factored expression is 

daab^rx / ,^2 



dp'rdy 



d^ab^[bb](Q)X' 

dQ'^dpj.dy 



^[bb]^rx{Q ) • (3) 
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FIG. 1: Hadronic production of an T via an intermediate heavy quark pair b and b. 

In writing Eq. Q, we approximate the transverse momentum and rapidity of the bb pair by 
Pt and y, respectively, because — 4ml ^ Vr- 

The function dd'ab^[bb]{Q)x' / dQ'^dp'^dy represents a partonic short- distance hard-part for 
inclusive production of a bb pair of invariant mass Q and quantum numbers [bb] . This short- 
distance hard-part is calculable in perturbation theory with the parton momenta of all light 
partons off-mass-shell by at least min(4m^,p^). The function J^[bb]^Tx{Q^) represents a 
transition probability distribution for a bb pair of invariant mass Q and quantum numbers 
[bb] to transmute into an T meson. It includes all dynamical bb interactions of momentum 
scale less than — 4m^. Different assumptions and choices for the transition probability 
distribution J-'{Q^) lead to different models of quarkonium production. We return to the 
topic of these models in Sec. IIIII 

The basic assumptions of this section imply that the transverse momentum distributions 
of the T states at transverse momenta pt ~ My will reflect the shape of the transverse 
momentum distribution for production of a bb pair whose mass Q ~ My. In this paper, we 
focus on the region below px ~ M^- If Pt ^ Q^, all final-state logarithmic terms of the 
form (a^ log(pf./(5^))^ can be resummed perturbatively to all orders in [2^. 

III. FIXED ORDER CALCULATION: pr ~ My 

When the transverse momentum, pj- ~ 0{My), the collinear factorized expression in 
Eqs. (j2I) and Q remains reliable with the partonic short-distance hard parts in Eq. (jS)) 
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computed as a power series in as in QCD perturbation theory. 

The transition probabihty distribution J^{Q'^) is introduced in Sec. |nj Different assump- 
tions and choices for J^{Q'^) correspond to different models of quarkonium production. In the 
color evaporation (or color-bleaching) model (CEM) 0, 21 1, an assumption is made, based 
qualitatively on semi-local duality, that one may safely ignore the details of the formation 
of color-neutral bound states with specific quantum numbers J^^ . In particular, in the case 
of states such as the J /if) and T that have J^^ = 1 , soft gluon effects are presumed to 
take care of whatever quantum numbers have to be arranged. Within our framework, this 
model is effectively represented by the statement that 

f Cr if 4m? < < AMI 
I otherwise . 

The non-perturbative constant Ct sets the overall normalization of the cross section. Its 
value cannot be predicted. It changes with the specific state of the T meson. In the CEM 
model, the parton-level T cross section in Eq. Q can be written as 

j-CEM Mil 



dp^dy JAml 



d(^ab^bb{Q) 

dQ'^dp'^dy 



(5) 



where the bb final-state includes a sum over all possible quantum states [bb] of the bb pair. 

In the color singlet model for quarkonium production 17|, a projection operator is used to 
place the bb system in the spin-state of the T, and explicit gluon radiation guarantees charge 
conjugation (C) and color conservation at the level of the hard-scattering amplitude. The 
distribution J^ibhi^riQ"^) is proportional to the square of the momentum-space wave function 
of the T, with the relative momentum of the bb pair = Q'^ — 4m^. Because the T 

wave function falls steeply, one can approximate ~ 4m^ in the bb partonic cross section. 
The integration / dQ"^ J^^j^i^_,-^{Q'^) in Eq. Q leads to the square of the T wave function at 
the origin |\I^(0)p. 

The non-relativistic QCD model (NRQCD) P, takes into consideration that the ve- 
locity of the heavy quark b (b) in the rest frame of the bb pair is much less than the speed 
of light. The velocity expansion translates into statements that the distribution J^(Q^) is a 
steeply falling function of the relative heavy quark momentum, = Q"^ — 4m^, and that its 
moments satisfy the inequalities 

{{qT) = / dQ' (g^)^^[,5HT(Q^) « {^rnir , (6) 



for moments, > 1. Correspondingly, one can expand the partonic hard part in Eq.Q at 
= (2mfe)^ and obtain 



7;iNRQCD 



dQ'^dp^dy 



= M^)J j dQ':F^,^^M') + O , (7) 



with mfc = Mt/2. The integral / dQ"^ J-'\p^^'^{Q'^) = (O[fc5](0)) corresponds to a local matrix 
element of the hh pair in the NRQCD model. In the NRQCD approach to heavy quarkonium 
production, the 66 pair need not have the quantum numbers of the T. It is assumed that 
non-perturbative soft gluons take care of the spin and color of the T. The sum in Eq. ((7|) 
runs over all spin and color states of the hh system. 

For the purpose of calculating the inclusive pt distributions of S-wave bound states 
at large enough py, both the CEM in Eq. and the leading order NRQCD approach 
in Eq. (|7j) are expected to yield distributions similar in shape because of the relatively 
weak dependence of the partonic hard-part in the limited range of . For example, 
Pt distributions of J/ip and ip' production at Tevatron energies are consistent with both 
CEM [2] and NRQCD calculations for > 5 GeV. For production of T{nS) states, we 
choose a J-'{Q'^) that covers both the CEM and main properties of the leading order NRQCD 



treatment of heavy quarkonium production 



2^ 



J' [f,b]^T(n5)lV ) — \ \°) 

I otherwise 

with z = {Q"^ — M^(^,^s))/{4:M^ — M^-^^^^). In Eq. (jH)), Cr(nS) and Q;x(n5) are parameters 
determined from data as discussed in Sec El With the choice of T^Q"^) in Eq. (jSJ, we 
reproduce the CEM by setting aT{n5) = and replacing the lower limit M'^^^g-^ by 4m^. 
Other than the color degree of freedom, we could mimic the features of NRQCD by choosing 
a very large value for aT(n5)- 

With our choice of J-'{Q'^), the transverse momentum distribution of T production be- 
comes 



dpj^dy 



Cr(nS) I ^ dQ 

J Ml, 



T(nS) 



AB^bb(Q)X 



(n5) 



4M| - M|(„^) 




dQ'^dpj'dy 

The hh cross section is factored in terms of parton densities and the partonic cross section 
as 

^^Q-I^ldy =^J^^^ M-^^ dx, M-,) . (10) 



The sum J2a,b runs over gluon and light quark flavors up to and including charm. The 
partonic cross sections, dcr^^_,f,ii^Q-^x / dQ^dp'^dy are computed at 0{al) from all 2-parton to 
3-parton Feynman diagrams for the subprocesses qq bbg, qg bbq, and gg — > bbg, with 



the squared amplitudes summed over the spins and colors of the bb pair 



24|. 



IV. THE REGION OF SMALL TRANSVERSE MOMENTUM 

When pt (or Qt of the bb pair) becomes small, the perturbatively calculated hard-part 
daah-^[bh\(Q)x' / dQ'^dp^dy in Eq. Q becomes singular 

dQ'^dp^dy p\ ' 

The I/pt singularity arises from the coUinear region of initial-state parton splitting. Gluon 
radiation from the final-state heavy quark lines does not contribute a collinear singu- 
larity because the heavy quark mass regulates this singularity. However, this gluon radiation 
does lead to a l/p% infrared divergence which should be absorbed into the non-local transi- 
tion probability distribution J^{Q'^) When p\ q^ = — Ami, soft gluon interactions 




FIG. 2: Diagram that illustrates multiple gluon radiation from an initial-state parton. 

between the spectator partons in the beam jets and the partons in JF(Q^) most likely break 
the factorization expressed in Eqs. (0) and 0- In this paper, our principal interest is to 
investigate how the large logarithmic terms from the initial-state gluon shower modify the 
I/Pt distribution when p'^ <^ M^. 
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A. Resummation of Sudakov logarithms in &-space 



Additional gluon radiation from the initial-state partons, recoiling against the bb pair as 
shown in Fig. |21 leads to (Sudakov) logarithmic contributions of the form log^(Q^/Pr) for 
each gluon radiation [2^. The effects of the large Sudakov logarithmic contributions, yery 
important in the region of small px, can be resummed to all orders in as when pr ^ Q j26l |. 
The resummation procedure tames the divergence seen in Eq. Adopting the Collins, 

Soper, and Sterman (CSS) impact-parameter 6-space (Fourier conjugate to pt) approach Q], 
we write the resummed transverse momentum distribution for bb production as 



db 



= J Jo{PTb)bWAB^[bb]{Q){b,Q,XA,XB) . (12) 

The function W^b^[{,5](q)(6, Q, xa, a;^) resums to all orders in QCD perturbation theory 
the singular terms from initial-state gluon showers that otherwise behave as 5'^{pt) and 
{l/p\) log™'((5^/PT) for all m > 0. In Eq. ()12|1 . the fractional partonic momenta are xa = 
and xb = with the overall center-of-mass collision energy, and y the 

rapidity of the hb pair; xa and Xb are independent of the transverse momentum px of the 
pair. The entire dependence on px appears in the argument of the Bessel function Jq. 

The expressions for the lowest order subprocesses gg —>■ bb and qq bb are independent 
of Pt- Therefore, the finite lowest order partonic cross sections can be used as prefactors in 
the overall 6-space distribution functions 2^. We write 

^AB^[bb]{Q){b,Q,XA,XB) = J2^gQi^^Q'^A,XB) ^^^'^''^^^^ 



dQ^ 
[b 

dQ' 



da^^"^^ - 

+ W-,,(6,Q,xa,xb)^5«. (13) 



The sum runs over all flavors of light quarks in the initial state. The lowest order partonic 
cross sections in Eq. (fT^. da^^^y^^{Q) / dQ^ with ij = qq,gg, depend on the choice of the 
production model. For our choice of J^bb^TiQ^)^ they are 



dQ' tiT dQ^ Q 



with 



£ -^"irw) (14) 

m 



ar(Q^)-'— 



9 3g2 



1 ■ 

1 + o7 



1-7; 



and 7 = 4m^/Q^. 

When the impact parameter b lies in the region much less than 1 GeV~^ where pertur- 
bation theory apphes, the distributions Wqq{b,Q,XA,XB) and Wgg{b,Q,XA,XB) in Eq. (fT!?jl 
can be expressed as j6| 

* (6, Q, XA, xb) = e-^-(^'«) f,/A (xA, ft, ^ f,/B {xB, fi, ^) , (16) 

where ij = qq and gg, fi is the factorization scale, and c = 2e~^'^ = 0{1), with Euler's 
constant 'Je ~ 0.577. All large Sudakov logarithmic terms from log(c^/6^) to log((5^) are 
resummed to all orders in as in the exponential factors with 



Sij{b,Q) 



In ( ^1 Aij{as{ii)) + Bij{as{ii)) 



(17) 



The functions Aij and Bij may be expanded in perturbative power series in a^; Aij = 
X]n=i -^Ij^ (c^s/^)") ^iid Bij = J2n=i^i^\'^s/'^)"'- The first two coefficients in the power 
series for Aij and the first term in the series for Bij are process-independent. For ij = qq they 
are the same as those that appear in resummation of the transverse momentum distribution 
fo. „,a.,ve-lepto„-pa. p..oduct,o„ (D..e„-Ya„ p.oduct,o„) fl Q. Fo. ^ they a.e 



the same as the coefficients that are appropriate for resummation of the pt distribution of 
H,ggs boson production HQQQ. 

The modified parton distribution functions in Eq. (fTH|) are expressed as |6[ 



In Eq. (IIHl), 4>a/A are the usual parton distribution functions. The functions Ca^i = 
J2n=oCli^i {oLsl'^Y 6-space coefficient functions with the lowest order terms normahzed 
to Cj^}^i{z, fi,c/b) = SaiS{l — z). All higher order coefficient functions are computed per- 
turbatively from the Fourier transform of the singular terms in p^-space from initial-state 
gluon showers with = c^/b"^. 

In the CSS formalism, Hij = 1 in Eq. All coefficient functions, C^^j with n > 1, are 
process- dependent Tepresenting the non-logarithmic short-distance partonic contributions to 
^resum_ Alternatively, one may be able to reorganize Eq. ()16|) such that all process- dependent 
short- distance contributions are moved into the process-dependent hard part Hij, leaving 



the coefficient functions Ca->i and the modified parton distributions process-independent . 
Expressions for if^" with n > 1 depend on the "resummation scheme", the choices made 
when the process- dependent finite pieces are moved from the higher order terms in the ,4's, 
i3's, and C's to the Hij functions jl2|. 

For production of a colorless object, such as the W , Z, and Higgs bosons or a virtual 
photon in the Drell-Yan process, all resummed {l/p^)\o^{Q'^ /p\) singular terms arise 
from initial-state gluon showers. For the T, which is not produced directly in the hard 
collision, additional singular l/p^^ terms can originate from soft gluon radiation from the 66 
pair. These additional l/pr terms should not be included in the resummation of Sudakov 
logarithms from initial-state gluon showers. The calculation of the n > 1 corrections to the 
coefficient functions C^^j and hard part h\^'' (or C^^^^ in the CSS formalism) should involve 
a systematic removal of these additional l/pr singular terms. 

The long-distance nature of soft-gluon radiation means that the additional singular terms 
from final-state radiation should be included in the non-perturbative transformation of the hh 
pair to the T. Therefore, the removal of the additional 1 /p\ terms depends on the models of 
T production. In terms of the two-step factorization procedure discussed in Sec. m the l/pr 
terms should be absorbed into the transition probability distributions J^[hb]^i:x defined in 
terms of matrix elements of non-local operators. In the NRQCD model of heavy quarkonium 
production, these l/p^ singularities are a consequence of the kinematic end-point of the 
quarkonium transverse momentum spectrum. Although the kinematic effect of soft-gluon 
emission from the heavy quark pair is usually a higher-order effect in the non-relativistic 
expansion, the high order non-perturbative contributions are enhanced in the region of the 
kinematic end point as 0, leading to a breakdown of the NRQCD expansion and the 

introduction of "shape functions" 27|, |28| . 

Faced with the model dependence of the corrections and the complication of separating 
singular terms of different origins, we resum only process-independent\ogajiihm\c terms from 
initial-state gluon showers in this paper. That is, we keep only A^^^g, -^g!]' ^^"^ '^gla 
Sudakov exponential functions Sq^g{b,Q) in Eq. (fTTjl . the lowest order coefficient function 
C^l^i, and the lowest order short- distance hard parts H-j^ = 1 in Eq. (fT^ . We choose the 
factorization scale fi = c/b for the resummed 6-space distribution in Eq. ()16|) 31[. We defer 
to a future study the calculation of process- and model- dependent higher order corrections 
C^^j and H^^^ as well as A^^^ and B^'^K As a result of these restrictions, we anticipate that 
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our calculation will somewhat underestimate the magnitude of the differential cross section 
in the region of small pt, and we return to this point in next section. 



B. Predictive power 

The predictive power of the Fourier transformed formalism in Eq. |T2|l depends critically 
on the shape of the 6-space distribution function bW{b, Q, xa, xb) 8]. Indeed, the resummed 
calculation of the transverse momentum distribution at low pt can be reliable only if the 
Fourier transformation in Eq. (fT^ is dominated by the region of small fe, where perturbation 
theory applies, and is not sensitive to the extrapolation to the region of large h. This 
condition is achieved if the distribution bWib.Q^XA.XB) has a prominent saddle point for 
6sp < 1 GeV-\ 

The location of the saddle point in the 6-space distribution depends not only on the value 
of Q but also strongly on the collision energy \fS (or, equivalently, on the values of the 
parton momentum fractions Xa and x^, that control the cross section) At 75 = 1.8 TeV, 
the saddle point can be within the perturbative region (fogp < 0.5 GeV^^) for Q as low as 
6 GeV 0. 

An extrapolation into the region of large h is needed in order for us to perform the Fourier 
transformation to the p^ distribution in Eq. (|12p. We choose the Qiu-Zhang prescription 
which has the desirable property that it separates cleanly the perturbative prediction at 
small h from non-perturbative contributions in the large h region. 

' ^ - (19) 

^ij^^ibynax, Q, XA, Xb) F^-^ib, Q] bmax) b > bmax 

for ij = qq and gg. The perturbative distribution Wlj^^^{b,Q,XA,XB) is given in Eq. (fTB|l . 
The nonperturbative function in the large b region is 



^max 



FS' = exp { - ln(^^) [g, {{b^ - (&LJ") + 92 {b' - bl 

-92 ffe' - O V (20) 



The (6^)" term with a < 1/2 represents a direct extrapolation of the resummed function 
Wfj^'^{b,Q,XA,XB)- The parameters, gi and a are fixed from Wf^"^^ if we require that the 
first and second derivatives of PFjj(6, Q, x^, x^) be continuous at 6 = b^ax- The xa and xb 
dependences of the nonperturbative function F^-^ are included in the parameters g\ and a. 
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The two terms proportional to 6^ correspond to power corrections in the evolution equa- 
tion. The g2 term represents a power correction from soft gluon showers. The g2 term is 
associated with the finite intrinsic transverse momentum of the incident partons. Since gi 
and a are fixed by the continuity of the Wij{b,Q,XA,XB) at 6 = bmax, the (in) sensitivity of 
the calculated px distribution to the numerical values of bmax, 92, and g2 is a good quanti- 
tative measure of the predictive power of the resummmation formalism j8|. 

In Fig. El we show the 6-space distributions that result from Eq. (fTI?|) at rapidity y = 
0. The functions are integrated over the mass range 2mb < Q < 2Mb- In evaluating 
the perturbative distribution, we keep only the process-independent terms in the Sudakov 
exponential functions Sg^g{b,Q) in Eq. (fTTj) and use the CTEQ6M parton densities 2^. For 
the extrapolation to the large b region, we choose bmax = 0.5 GeV~^ and g2 = g2 = 0. The 
magnitude of bWgg{b, Q, xa, xb) is scaled by a factor of 20, as is indicated in the figure. The 
gg contribution far exceeds the qq contribution to T production at Tevatron energies. For 
both channels, the saddle points are clearly defined and have numerical values well within 
the perturbative region. For the dominant gg channel, the location of the saddle point is 
bsr, ~ 0.25 GeV^^. For perspective, we remark that this value is smaller than that for the 
saddle point of W and Z boson production at Tevatron energies |8|. This feature arises 
because the gluon anomalous dimension at small x is much larger than that of the quarks, 
compensating for the fact that the mass Q here is much less than the mass of W and Z 
bosons. This analysis leads us to expect that the QCD resummed pt distribution of T 
production in Eq. ()12|1 can be predicted reliably in the region of small and intermediate pt 
because it is dominated by perturbative contributions in the region of small b. 



V. NUMERICAL RESULTS 



In this section we present the results of our numerical computation, including a compari- 
son with data. We compute the T transverse momentum distribution from Eq. Q. For the 
region of large transverse momentum, px ~ 0{Mx), the bb cross section is given in Eq. ()10p. 
For the region of small px, we use the all-orders resummed bb cross section in Eq. p2|) with 
the 6-space distribution WAB^bb specified in Eqs. f|T3|l and (fT9|l . We set rjif, = 4.5 GeV, 
and we use a two-loop expression for a^, in keeping with our use of the CTEQ6M parton 
densities 12911. 
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0.5 1 1.5 2 0.5 1 1.5 2 

b{1/GeV) b{1/GeV) 



FIG. 3: The 6-space distributions for T production: (a) gg channel, and (b) the sum of all qq 
channels. Note that the magnitude of the gg distribution has been scaled by a factor of 20. The 
functions are evaluated at rapidity y = and are integrated over the mass range 2mf, < Q < 2Mb- 

To distinguish the production of T{nS) states with different n, we choose different powers 
<^r{nS) and normalization constants Cx(n5); in addition to the differences in mass threshold 
on the limits of the dQ^ integration in Eq. Q. The values of ar(nS) and Cr(nS) are correlated. 
A larger value of ax(n5) leads to a larger value of Cr{nS)- 

A. Matching of results at small and large pT 

In a complete calculation, one would expect a seamless joining of the results applicable at 
small and at large pt- In the CSS resummation formalism for production of a color singlet 
heavy boson, this matching is accomplished through the introduction of an "asymptotic" 
term, a^^^, and 

da = da™' + {da^"'' - da'^y'^) . (21) 

The term cr'^'^^™ is constructed to cancel the singular behavior of a^'^''* as px ^ and to 
cancel cr^*^'^"™ when pt ^ Q- It is obtained from the fixed-order terms in the expansion of 
^resum ^ power scrics in a^. 

The procedure just described is not immediately applicable in our case. Because the bb 
system is not necessarily in a color neutral state, ctP'^''* includes radiation from the heavy 
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quark system as well as from the incoming partons. This final state radiation is not included 
in either a^'^^^^ or cr'^'^^™. Soft gluon radiation from the heavy quark system leads to an 
infrared divergent l/px singularity [l^ which should be absorbed into J^{Q'^) To avoid 
extrapolation of ctP'"'* into region of low p-p, we adopt the following matching procedure: 

{J rcsum 
4B->T(n5)X 
d^p^ Vt < Vtm 

Matching is done at a value pt^ chosen as the location of intersection of the resummed 
and perturbative c omp onents of the Pt distribution. From other work on resummed Pt 
spectra , we expect ~ Mr/2. 

To ensure a smooth parameter-free matching of cr^^^'^™ to the perturbative pt distribution 
^pert computed at 0{al), we would need to calculate the process-dependent O (as) corrections 
C^l,i and H^j^ (or C^^j in the CSS formalism) for 0-'''='^"™. If these 0{as) corrections were 
included, a'^^^^'^ would also be of order 0{a^) at the matching point where the logarithms 
are not important. Based on prior experience |1^J, we expect that these effects will increase 
the predicted normalization of da^'^^^^ / dp\dy , and change the shape of the px distribution 
somewhat, increasing (decreasing) the spectrum at small (large) px- For the reasons stated 
in last section, we do not calculate the order corrections to 0"'''='^"™ in this paper. To 
account for the size of the order corrections, we introduce a resummation enhancement 
factor Kr such that 



■^(0) ^(1) Os' 



MO) ,Mi) 



VT . 



^.Cll^Cll^ffif. (23) 



We assume that Kr is a constant. The factor Kj. should not be confused with a "i^-factor" 
for the overall px distribution. It is invoked because we do not calculate the order as 
corrections to a"^^^^™, and our ci'^^'^"™ ~ C^(«s) when px ~ 0{Mr)- 

Displayed in Fig. |3 are curves that show the differential cross section for production 
of the T(IS') as a function of px- The curves in the region of large px illustrate the de- 
pendence of the fixed-order 0(«^) perturbative cross section on the common renormaliza- 



tion/factorization scale /i. We vary /i over the range 0.5 < /i//io < 2 where /xq = yQ"^ + Px- 
This variation demonstrates the inevitable theoretical uncertainty of a fixed-order calcula- 
tion. It could be reduced if a formidable 0{a^) calculation is done in perturbation theory. 



We fix /X = O.SyQ^ + Px remainder of our discussion. The 1/px divergence men- 

tioned in Eq.()ll|) is evident in the fixed-order curves. Shown for purposes of perspective is 
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a dot-dashed line that represents the resummed prediction, apphcable at small pt, obtained 
with bmax = 0.5 GeV~^ and g2 = 92 = 0. The set of curves illustrates the range of possibil- 
ities for the value of the matching point Ptm- To obtain these results, we set Cx = 0.044, 
ar = 0, and Kj. = 1.22, for reasons that are explained in the next subsection. 




PT(GeV) 



FIG. 4: Inclusive transverse momentum distribution of the T(IS') in pp interactions at ^/S = 
1.8 TeV. The solid, dashed, and dotted lines in the region of large pT are obtained from fixed-order 
0{al) perturbative QCD for three different values of the scale solid for fj, = 0.5/xo, dashed for 
fi = fiQ, and dotted for p = 2fiQ; with fiQ = \J Q"^ + Pt- The dot-dashed line in the region of small 
Pt is the resummed prediction. 



B. Comparison with data 

In order to make contact with data we must determine values for Kr in Eq. (j23p and 
PTa/ in Eq- (1221); and for the two non-perturbative parameters Cr and ax in the transition 
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probability distribution JF(Q^). The structure of Eq. indicates that the data at highpT 
determine Cr and ar, and the data at lowpx fix the enhancement factor Kr- Dependence on 
the parameter a-^ turns out to be very weak, as might be expected from the hmited range 
in over which J-'{Q'^) is probed. This weak dependence confirms that the production 
models we consider predict very similar inclusive pt distributions. We choose to set ar = 
for all three T states. Second, common values of Kr and pt^^^ work adequately for all three T 
states, as might be expected since the differences are small among the three Upsilon masses. 

In our approach to the data, Cx represents not just the normalization in the transition 
probability distribution JF(Q^), but the product of this normalization times the unknown K- 
factor from order 0{a'^) perturbative contributions at large pt- It should not be surprising 
that non-perturbative free parameters enter the comparison with data at large pt- The 
reliance on non-perturbative parameters to set the normalization is true of all models other 
than the color-singlet model 

We compare our calculation with data published by the Collider Detector at Fermilab 
(CDF) collaboration |4,Q] obtained in run I of the Tevatron collider at y/S = 1.8 TeV. In 
the second of the two publications, it is noted that the measured rates are lower than those 
reported in the first paper. To account for the difference in our fits to the data, we include 
an overall multiplicative normalization factor C„, whose value we determine from our 
fitting routine. This factor is used only for the 1995 data P]. 

Following our initial qualitative exploration of the data, we are left with the three pa- 
rameters Cr, the common values of Kr and Ptm^ data adjustment factor We 
use a minimization procedure to determine these quantities. We find best fit values 
Kr = 1.22 ± 0.02, and p^Af ~ 4.27 GeV. The value of Kr is comparable to typical K-factors 
found in next-to-leading order calculations, but, as remarked above, the origin here is dif- 
ferent. The matching point is fixed essentially by the location where the resummed and 
perturbative cross section intersect. Its val ue, vr.. ~ Mr/2, is similar to results found in 
other work on resummed px spectra ll|, 112, |13| . We find that the values of Cr- 
0.044 ± 0.001, 0.040 ± 0.006, and 0.041 ± 0.003 for T{1S), T{2S), and T(3S), respectively, 
are approximately independent of Mr, meaning that the differences in rates for the three 
5- wave T states are accounted for by the different threshold values of the integrals in Eq. Q. 

In our determination of Cr, Kr and Ptm^ we keep only process-independent terms in 
the Sudakov exponential functions Sij{b,Q) in Eq. (fTTj). and the parameters of the non- 
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perturbative function F^^ are fixed at bmax = 0.5 GeV~^, and g2 = 92 = 0. Because of 
the dominance of the perturbative small-6 region under the curves of bWij in Fig. El any 
reasonable values of g2 and g2 lead to transverse momentum distributions that do not differ 
more than one percent from those calculated with g2 = 92 = ^ l^^^- Without adjusting the 
normalization, we find a few percent change in the resummed distributions over the entire 
low Pt region when we vary bmax from 0.3 to 0.7 GeV~^. 

The principal predictive power of our calculation is the shape of the pr-distribution for the 
full Pt region. In Fig. we present our calculation of the transverse momentum distribution 
for hadronic production of T{nS), n = 1 — 3, as obtained from our Eq. (j221), and multiplied 
by the leptonic branching fractions B. We use the values of B from Ref. IsJ]. The solid 
lines are for b^ax = 0.5 GeV^^ while the dashed and dotted lines are for bmax = 0.3 and 
0.7 GeV~^, respectively. Also shown in Fig. El are data from the CDF collaboration 
We determine a data normalization adjustment of C„ = 0.88 ± 0.05 and use this value to 
multiply only the 1995 cross sections shown in the figure. The shapes of the pt distributions 
are consistent with experimental results. 




FIG. 5: Calculated differential cross sections times leptonic branching fractions B, evaluated at 
y = 0, as functions of transverse momentum for hadronic production of (a) T{1S), (b) T(2S), and 
(c) T(35), along with CDF data P, ^ at = 1.8 TeV. The dashed, solid, and dotted lines show 
the result of our full Cctlculs-tion for huiax — 

0.3, 0.5, and 0.7 GeV-^ respectively. The 1995 CDF 

cross sections are multiplied by C„ = 0.88. 

The essential similarity of the production differential cross sections for the three T{nS) 
states is illustrated in Fig. IHl Shown are the differential cross sections divided by their re- 
spective integrals over the range < < 20 GeV. The integrated values are computed from 
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the theoretical cross sections and used to scale the experimental as well as the theoretical 
results. The three theory curves are practically indistinguishable. The transverse momen- 
tum distribution is described well over the full range of pt- Since the curves in Fig. IHl are 
normalized by the integrated cross sections, dependence on the normalization parameters 
Ct cancels in the ratio. The shape for < My/2 is predicted quantitatively. It reflects 
the resummation of the gluon shower and is independent of parameter choices. The good 
agreement with data over the full range in pt is based on the choice of only two adjustable 
constants, the resummation enhancement factor Kr and the matching point Ptm- 
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FIG. 6: Normalized transverse momentum distributions for T production: T(IS') (solid), T{2S) 
(dashed), and T(35') (dotted), along with the 2002 CDF data [ij. 



18 



VI. SUMMARY AND CONCLUSIONS 



In this paper we calculate the transverse momentum pt distribution for production of 
the T states in hadronic reactions, applicable over the full range of values of pt- Our 
starting assumption is that the pt distribution of T production may be derived from the pt 
distribution for the production of a pair bb of bottom quarks. We express the differential 
cross section in terms of a two-step factorization procedure. We justify the validity of an 
all-orders soft-gluon resummation approach to compute the px distribution in the region 
of small and intermediate pt where pt < Mx. Resummation is necessary to deal with the 
perturbative l/p^ singularity and the large logarithmic enhancements that arise from initial- 
state gluon showers. We demonstrate that the px distribution at low px is dominated by the 
region of small impact parameter b and that it may be computed reliably in perturbation 
theory. We express the cross section at large px by the 0{a^) lowest-order non-vanishing 
perturbative contribution. Our results are in good agreement with data from pp interactions 
at the Fermilab Tevatron collider at center-of-mass energy \^ =1.8 TeV, and they confirm 
that the resummable part of the initial-state gluon showers provides the correct shape of 
the Pt distribution in the region of small pt- 

An improvement of our calculation in the region px < Mr would require inclusion of the 
order as process-dependent corrections associated with the coefficient functions C^^-* in the 
CSS formalism (or equivalently, C'^^^ and H-j'' in Eq. p6|)). Based on prior experience jisl ]. 
we expect that these effects will increase the predicted normalization of da^^^™ / dpTdy , and 
change the shape of the pt distribution somewhat, increasing (decreasing) the spectrum 
at small (large) px- A complete calculation of the order corrections C^^j and iff^^^ in 
Eq. (fTBj) would provide a better test of QCD predictions. In the region of large pr, an 
improved prediction of the normalization and shape of the differential cross section would 
require a formidable O(a^) calculation of bb production. 

Inclusive production of the T states in the central region of rapidity at Tevatron energies 
and above is controlled by partonic subprocesses initiated by gluons. The typical value of 
the fractional momentum x carried by the gluons is determined by the ratio Mx / V^- The 
growth of the gluon density as x decreases leads to two expected changes in our predictions 
for larger v^- First, and perhaps obvious, the magnitude of the cross section near the 
peak in, e.g., Fig. El will increase. Second, and more subtle is the prediction that the peak 
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location should shift to a greater value of px as ^/S grows. We use the same parameters 
as those at ^/S = 1.8 TeV, Fig. El to compute T(15') production at the Tevatron in run-II 
at "s/S* = 1.96 TeV; our results are shown in Fig. [7| The change of y/S from 1.8 TeV to 
1.96 TeV does not produce a marked difference in the spectrum, but we expect the shift of 
the peak in px to be about 1 GeV at the LHC energy of y/S = 14 TeV. 




FIG. 7: Differential cross sections times leptonic branching fractions B, evaluated at y = 0, as 
functions of transverse momentum for hadronic production of the T(IS') at "v/S* = 1.8 TeV (solid 
line) and 1.96 TeV (dashed line). 

Our focus on T production may motivate questions about the analogous production 
of the J/ip states. The mass of the J/tp is relatively small, meaning that inverse power 
contributions of the form 1/Q are potentially as significant as the logarithmic terms log{Q) 
that we resum. In addition, the saddle point in the 6-space distribution moves into, or close 
to the region in which perturbation theory can no longer be claimed to dominate the pt 
distribution. 
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